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We consider metallic carbon nanotubes with an overlying unidirectional electrical chiral (wavevec- 
tor out of the radial direction, where the axial direction is included) superlattice potential. We show 
that for superlattices with a wavevector close to the axial direction, the electron velocity assumes 
the same value as for nanotubes without superlattice. Due to an increased number of phonons with 
different momenta but lower electron-phonon scattering probabilities, we obtain a large enhance- 
ment of the high-voltage conductance and current sustainability in comparison with the nanotube 
without superlattice. 

PACS numbers: 63.22.Gh, 73.21. Cd, 73.63.Fg, 73.50.Fq 



Depending on their chirality, carbon nanotubes (NT) 
behave either like a semi-conductor or a metal. In the 
first case, they offer interesting alternative for building 
logical circuits. In the second case, they can be used as 
nanometer-sized metallic wires in logical circuits. This is 
particularly useful since they can sustain very high cur- 
rents before breaking. At low voltages (U < 0.17V) the 
effective electron scattering length at room temperature 
in metallic NTs is mainly governed by acoustical phoiion 
and impurity scattering with a value of a few hundred 
nanometers [l|. At higher voltages, scattering with hot 
optical phonons created by electron-phonon scattering 
becomes relevant. This leads to a significant reduction of 
the electron's mean free path down to roughly l sc w 10 
nm resulting in a large increase in the absolute 

and differential resistance. Due to the large number of 
optical phonons, phonon-phonon scattering with acoustic 
phonons produces heat in the NT that ultimately causes 
the electrical breakdown 0, [I[ . 

In Ref . it was argued that the performance of a metal- 
lic NT, i.e., its absolute and differential conductance, can 
be enhanced considerably by isotopical disorder enrich- 
ment. This causes additional relaxation paths for opti- 
cal phonons by disorder scattering. The purpose of the 
present letter is to propose a different mechanism to 
enhance the electronic transport. We show that by ap- 
plying an unidirectional electrical superlattice (SL) (cf. 
Fig. 1) with wavevector close to the axial direction of the 
NT, we can enhance the (differential) conductance con- 
siderably, especially in the large voltage regime. Such a 
potential could be for example produced by adatom de- 
position via electron beams directed on the NT [l(| or (at 
least approximately) by twisted periodical patterned top 
and bottom gate electrodes or the coupling of the NT to 
surface acoustic waves [Til ]. Since the (average) phonon 
number is proportional to the inverse electron-phonon 
scattering time l/r op , and the inverse electron mean-free 
path l/l sc is proportional to the phonon number times 
1 /r cp in the hot phonon regime, we obtain a quadratic de- 
pendence of the electron mean-free path on the scattering 



time ^ sc ~ Tg p . Below it will be shown that an applica- 
tion of an electrical chiral potential causes a large number 
of different phonons to take part in the electron-phonon 
scattering process with increased electron-phonon scat- 
tering times, so that l/r cp ~ 2jV T ep- This is what 
causes the strong decrease of the (differential) resistance 
that scales with l/l sc ~ J2i^/( T c P ) 2 ^ V T o P by using 
Matthicssen's rule and reduces the phonon temperature. 

It was shown recently that new Dirac points in the 
energy spectrum can be opened by imposing an SL on 
the graphene lattice [L2h16| ] . This is also seen in NTs for 
potentials with wavevector in the radial direction. Wc 
will show that they vanish for general chiral potentials. 

The Hamiltonian near the Dirac point K for a NT with 
axis in y' direction subjected to an SL potential reads 
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V(x> + t lV ') 
—ihvF {d x i + idyi) 



-ihvF (d x i — idyi) 
V(x' + t y y') 
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where vp is the Fermi velocity. For a NT with circumfer- 
ence D, the SL potential V is periodic in the radial (axial) 
direction with periodicity d (d/t 7 ), i.e., V(x'+d) = V(x') 
and one has D/d € N. In the following we solve the 
eigenvalue equation Hku(t') = eu(r'). We use the ab- 
breviation try = tan(7) where 7 is the chiral angle of the 
SL potential V. The metallic NT boundary conditions 
are given by u(x' + D, y') = u(x' , y') [l~7l ]. 

To solve the eigenvalue equation we follow first a 
transfer matrix method similar to Ref. UH By us- 
ing the coordinates x = x' + t 7 y' and y = y', the 
solution of the Schrodingcr equation has the Bloch 
form u(r) = e lqv {u 1 {x) 1 U2(x)) T with (u 1 (x),u 2 {x)) T = 
A(x)(ui(0),u 2 (0)) T , where 



A(x) = e- iqx T exp 

q/Tr 



dx'M 



V (a 
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V(x) 



iK(x)/(l + ity) . 
iK{x)/{\ - ity) -q/T* , 
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and q — qty/Tz, k(x) = [e — V(x)]/Hvf, 
The operator V indicates path ordering and places all 



t 2 . 

7 
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larger values of x to the left. The Bloch condition 
reads (ui(d), u 2 (d)) T = e l,? (ui(0), u 2 (0)) T with det[e i?? - 
A(d)] = when d is the SL wavelength. 

Consider first a carbon NT in a chiral periodical lattice 
of two piecewise constant potentials of the form 



V(x) 



Vi 

v 2 



if < x < di 
if g?i < x < d\ + d 2 . 



(4) 



Then we obtain A(c?) = A1A2, where d = d\ + d 2 

A, : = e 



Iccsla.m + ^^MyM, (5) 



Ot e (di) 



a e [d i ) = d iy li$/T*-<P/T*. 



(6) 



Ki is given by k(x), with V(x) = V*. Since dct[c lqx A(x)] = 
1 we have for the eigenvalues of the matrix A(d), £ = 
e- i ? d (l/2)(T± VT 2 -4), with T = Tr[e i£ i d A(d)] and 

T = 2 cos[a(cq)] cos[a(d 2 )] (7) 



sin[a(di)] sin[a(d2)] / q 2 n\n 2 



a(di)a(d 2 ) 



T 4 
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By taking into account that T is real, we obtain from ([7]) 
the dispersion relation 



2 cos{qd + 77) = T . 



(8) 



In the following, we restrict ourselves to mirror sym- 
metric potentials V\ = — V 2 = V, with d\ = d 2 = d/2 
leading to the best current-voltage results over all two- 
step potentials. This leads to an energy spectrum that 
possesses a mirror symmetry at e = as a function of the 
quasi-momentum q in y-dircction. For e = 0, we obtain 
from Eq. ® that T = 2 + q 2 d 2 sm 2 [a (d/2)}/al(d/2)T^. 
This leads with (J5J) to the existence of new Dirac points 
flij at zero chirality t 7 = 0. The number of these points 
is given by [V d / '2irhv fT 7 ] where [x] is the largest inte- 
ger number smaller than x. For t 7 7^ an energy gap is 
opened and the Dirac points disappear (see left panel in 
Fig. 1). 

Next we calculate the energy values of the bands at 
zero momentum q = 0. Eqs. ([7]) and ([8|) delivers for 
these energy values e''(0) = (±77 + 2im)T 1 hvF/d where 
n determines the energy bands for fixed quasimomentum 
77. Thus the energy bands are far more separated in en- 
ergy space for t 1 3> 1 than for the system without chiral 
potential, i.e. e v (0) for i 7 = 0. 

In order to see how the lowest energy band scales with 
V and t 7 , we calculate from ([7|) the energy dispersion 
of the lowest band for metallic NTs, i.e. 7/ = 0, in the 
regime |e s | < hv F T 7 /d, V and q 2 < T 2 (V(x)/hv F ) 2 , to 
be called 1Z. We then obtain 



e s =sriVF\ I h^i — h i—dr sin 
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FIG. 1: (Color online) Upper panel shows a NT with an 
overlying chiral electrical superlattice potential. Lower panel 
shows the two lowest energy bands (d = D) by solving Eq. (|8} 
for 77 = as a function of the rescaled axial quasi-momentum 
qd/2 for various chiral angles 7 with Vd/2hvF = 5 (left panel) 
and various chiral potentials Vd/2hvF with ty = 1 (right 
panel). Note that the dotted curve in the right panel crosses 
the x-axis exactly only at q = 0. 



where s = ±1 and sgn[x] is the sign of x. Note that 
T = sin[a (d/2)]e iQo(d/2) /a {d/2) for the symmetric two- 
step potential (|4|). 

In Fig. 1, we plot the two lowest energy bands (d = D) 
for 77 = 0, by solving Eq. (JSJ numerically. Eq. © leads 
to the electron velocity v y '(q) = de s /d(hq) along the NT 
axis (the q dependency of q has to be considered in the 
derivate). We restrict our discussion to momentum re- 
gion near the Dirac point, i.e., qd/2 < 1. It is smaller for 
larger potentials V being maximal at evo — > with value 
v y >(q) < vf- We point out that in general for |i 7 | > 1 
we have v y >(q) ~ vf in 1Z irrespective of the potential 
strength V. Since electron-phonon scattering times are 
proportional to v y i (q) we restrict our transport calcula- 
tions below to NTs with |£ 7 | > 1 leading to the largest 
conductivities. For the group velocity of the electrons in 
radial direction we have v x >(q) = de s /d(hq) — tyV y '(q). 
For t 7 > 1 we obtain v x ,(q) fav F (l- |r| 2 )/T 7 in TZ lead- 
ing to the collimination of the electron beam [lU, i.e., 
|wj,'(g)| > This expression is even valid where 

t 2 > |r| 2 . On the other hand for t 2 < |T| 2 we obtain 
v X ' w tyVp(l — |r| 2 )/r leading to the vanishing of collim- 
ination at chiral angles |r| 2 /(l - |r| 2 ) < |t T | < (1 - |r| 2 ). 
For even smaller i 7 we obtain collimination again. 

In Fig. 1 we see how the energy bands oscillate, in 
accordance with Eq. ([9]) for q < T 7 V/hvF, thus forming 
band SL valleys. The central valley possesses a true Dirac 
point at q = 0. The SL side valleys have then a minimum 
at sin 2 (<7<i/2) = 0. Within the SL valleys, electrons travel 
either to the left (right) for de s /dq < (de s /dq > 0). 
The number of side valleys can be read of from §§§ as 
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2m i with mi — [t 7 Vd/T 7 Hvf2tt]. 

For SLs with wavevector in exact axial direction we 
can read off the physics from the chiral case by choosing 



n 



where d is now the wavelength of 



the SL potential in axial direction with quasimomentum 
k. The wavevector q in the circumferential direction is 
quantized by q — lixnjD due to the periodic boundary 
conditions of the wave functions. We point out that also 
the wavcfunctions Eqs. (fTTj) . (fT2|) . and the considerations 
below Eqs. (fTTj) . (fTB| are still valid with the additional 
replacements k x — > k y and x — > y. Eq. (|9|) delivers that 
the energy bands are separated by HvFqT which means 
that the energy spacing between the energy bands goes to 
zero for infinite potential strength. From Eq. ([9|) we have 
v — > vf for Vd/hvp — ► oo where for the lowest band, i.e. 
q = 0, v — > vf even for finite potentials Vd/HvF- 

Next, we determine the eigenvectors v v of the matrix 
A(d). These are given by v r ' = ^-([a+i sm(r]+qd)]/b, 1) T 
with 



" j I \ sin ( a 2) , 7 / xSin(ai) 
qd2 cos(ai) h qd\ cos(a2) 

.sin(a 2 ) , .sin(ai) 

K2«2 cos(ai) YK\d\ cos(a2j 



1 

rp2' 

1 



a 2 



Ql 



. sin(ai) sin(a2) d\d 2 q (ki — k 2 ) 
aia 2 T^(l — it 7 ) 



1 — if- 



(11) 



and N\ is a normalization factor. The eigen- 
function u v (x, q) of the full Hamiltonian Hk 
is then given by u v (x,q) = (cos[ao(x)]/ + 
{sin[a (a:)]/ao(a;)}My)v'' for x < d/2. Eq. (fTTj) 
leads to (u''(a;, q)\e lk " :X I\u~ ,1 (x, — q)) = where we 
used the abbreviation (u q {x^q)\e lk ' cX a\u~ ri {x,—q)) = 

J dx (u , '(x,q)\e' lkx <7\u~ , '(x,-q)) with a G {cr x ,cr y ,/}. 
Here / is the identity matrix and a x , <r y the spin 
matrices. The wavevector of the phonons or impurities is 
denoted by k x . This is a generalization of the results that 
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inner- valley backward impurity scattering in Refs. [l9l 
and deformation potential phonon scattering in Rcf. 
does not exist in the lowest band in metallic NTs in 
contrast to semiconducting ones. 

We then obtain from for the lowest band eigen- 
functions u s (x,q) for rj = in the regime 1Z correspond- 
ing to the eigenvalues © 



u s ( X ,q) = l-e-^ 



x [ — ^ sm( qd) 
\qd 



./l - if 7 \ T, 



hvpq 



v/i 

<£*(x)- 



it- 



(12) 



i + 



zt- 



it- 



</>(x) 



where iY u in f| 12[) denotes a normalization fac- 
tor. The phase factor </>(x) is given by = 
exp [i J dx'sgn [^(x^JaoCx')/^'] . We point out that 
the lowest band eigenvalues ([9]) and eigenfunctions 
(|12p are more generally valid for chiral poten- 
tials V(x) = V(x + D), where we have to as- 
sume that /„ D dx'sgn [V(x')]ao(x')/x' = 0. In or- 



der to derive the eigenfunctions (fT2|) we first for- 
mulate the eigenvalue problem corresponding to (fTJ) 
in the basis (^1 - if 7 , yj\ + it^ e~' lqx <j)* {x)e iqv and 
{-yjl - it 7 , + it 7 ) 7 V*"? 1 '(j){x)e lqv '. The resulting 
Hamiltonian is evaluated perturbatively in lowest or- 
der in (hvF/d)^ sm(qd)a z and (?iUi?g/T 7 )(Re[0 2 (x)](Tj,-|- 
Ini[(/> 2 (x)]£r x ) resulting in © and (fT2")) when the ~ q 2 



terms in ao are absent. The full expressions ^ and (JT2 
valid in the regime 1Z can then be read off by comparing 
the first order expressions with the formal solution of (jSJ) 
for general chiral potentials V(x) leading effectively to 
the ~ q 2 correction factor in ao ©. 

By using (fT2|) we are now able to calculate the back- 
ward squared transition matrix elements in the regime 7Z 
being inverse proportional to the inverse electron phonon 
scattering time. Here we restrict us to SL inner and 
inter-valley backward scattering, i.e. d q sin(gd/2) 2 ^ 
and d q ' sm(q'd/2) 2 «s 0, at large chiral angle |t 7 | > 
1, for (qd) 2 T 2 < 4T 7 4 sin(gd/2) 2 and (q'd) 2 T 2 < 
4T 4 sin(<fd/2) 2 which is the relevant regime for transport 
at high applied bias voltages. We obtain 

\(u s (x,q)\ae^ x \u s ,(x,q'))\ 2 « (A° a )%^ (13) 
1 



2to- 



±2ir[2Vjd/2iThv F T,,] ' 



3=1 ± 



0. 



where (A°) 2 « and (^J 2 a 1, (^) 2 » (yL}) 2 

Here we introduced the abbreviation k x = k x — q' + 
q. Note that we used in (fTB")) the approximation that 
otoidi) ~ di\Vi\/hvp T-y valid for the majority of SL val- 
leys. 

For the symmetric two-step potential ((4|) we have 
7712 = 1 and V\ = V in (fTB)) . We generalized in (fTB)) 
our results to a chain of symmetric two-step potentials 
with d\ = d 2 = d/2ni2 of potential heights Vj where 
we assume the potential heights are separated consider- 
ably [2V l d/2Trhv F T 1 ] ^ [2V j d/2irhv F T 1 ] ^ for i j. 
This restriction implies that the number of different 
phonons taking part in an elcctron-phonon scattering 
process is maximal. Eq. (fTB)) delivers that for every 
phonon type of certain momentum the scattering prob- 
ability is l/2ui2 smaller than in the case of no existing 
chiral potential 0. By using (fT2)) we obtain for gen- 
eral step potentials within the same appromation used 
in (ED E fcx \(u s (x,q)\ae lk * x \u s ,{x,q'))\ 2 w 1 for a = a x 
and zero for a G {a y ,I}. This shows with the help of the 
discussion in the second paragraph, that the considered 
chain of symmetric two-step potentials (|13() should give 
the best transport results over all step potentials and that 
even a general step potential should show an enhanced 
conductivity compared to the pristine NT. 

The number of SL valleys for a chain of two-step po- 
tentials is now mi = [|t 7 |min[Vi]ti/r 7 fit;f 27r] which can 
be read off from Eqs. (f5)) and (0 as in the case of the 
symmetric two-step potential. Here we denote min[Vi] 
as the minimum of all Vis. Finally, we mention that for 
the forward scattering amplitudes, i.e. d q s'm(qd/2) 2 ^ 
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FIG. 2: (Color online) Upper panel: Current-voltage charac- 
teristic of the (mi, 77i2 ) model (2m i SL side- valleys and 2m2 
phonon species) for mi = 1 (green dotted curves), mi = 2 
(blue dashed curves) and mi = 3 (red dashed-dotted curves). 
Curves of the same style have different 7712 parameters with 
(mi, 1) (bottom curves), (mi, mi + 1) (middle curves) and 
(mi, (mi + l) 2 ) (top curves). This ordering is reversed in the 
inset and the lower panel. The NT length is L = 300nm. 
The black solid curves show the current-voltage characteris- 
tic with zero chiral potential V — and hot phonons (n* 7^ 0) 
or frozen phonons fixed at zero temperature (n 8 = 0) [f|. In- 
set: Differential resistance as a function of the NT length L. 
Lower panel: Position and energy averaged phonon number 
of SL inner- valley T- phonons n r for L — 300nm [g|. 



and <9 g < sin(g'd/2) 2 § 0, wc have (A°J 2 = (A") 2 = 1, 
(A° J 2 = and (Al ) 2 = in 03]). 

Until now, we have ignored transitions of electrons be- 
tween the K and K' valleys. For the eigenvalue problem 
of the K' valley we can repeat the discussion above for the 



K valley by using the substitution t 



-£ 7 and q 



in the corresponding expressions |2l( . Zone boundary A[ 
phonon backward scattering is the only relevant phonon- 
scattering mechanism in this case 0] . We now have 
to calculate the square of transition matrix element (|13|) 
with a = <jy [22I ] where one of the wavefunctions stands 
for a K' valley function and the other is an eigenfunc- 
tion of the K valley Eq. (13| leads to (A° y ) 2 « and 
{A^) 2 rj 1. The corresponding forward scattering am- 
plitudes vanish. 

Let us now consider the conductance of a NT with an 
SL potential. To this end we use the electron-phonon 
Boltzmann approach with the parameters established in 
Refs.lol.l6l.l23l for a NT lying on a substrate without chiral 
potential where also a thorough discussion on the scat- 



tering mechanisms can be found. The leads and the sub- 
strate is fixed at room temperature. Before going into the 
details of our calculation we want to recall here that the 
number 2mi was defined as the number of SL side- valley 
and 2?7i2 are the number of steps in the SL chain po- 
tential which is equal to the number of different phonon 
species taking part in a backward scattering process. 

In the following, we assume that higher bands do 
not contribute to the conductivity so that we have 
HvfT^it I D > eU (hvF2ir/d > eU) for wavevectors of 
the SL out of (exact in) axial direction. For large enough 
applied bias voltages eU S> |e s (<7i)| where qi is defined by 
sin(^d/2) = 0, q 2 < T 2 (min[ Vj]d/hvp ) 2 and large chi- 
rality, i.e. i 7 > 1, we obtain the following idealized band 
system: We have one central SL band with dispersion 
e(k) rj ±fi,Wi?|/c| and 2mi SL sidebands with momen- 
tum shifted dispersion e(k) rj ± lim^ _>o ^>f \Jk 2 + k 2 . 
Phonons of 2m2 types contribute to electron-phonon 
backward scattering with scattering times 2m2T^ p . Here 
v stands for T for longitudinal E2 zone-center scatter- 
ing or K for A[ zone boundary phonon scattering 0, [|| . 

is the corresponding electron-phonon scattering time 
without chiral potential. We can simplifyour calculation 
by using the same phonon velocity v£ p p, @ of the sys- 
tem without SL for all type of optical phonons. This is 
justified by the fact that our results do not depend much 
on the specific velocity value since the phonon mean free 
path is much smaller than the NT length as we have ver- 
ified numerically. 

It is enough to consider only forward scattering be- 
tween the central SL valley and the 2mi side valleys 
mediated by transversal optical T phonons with scatter- 
ing time rjp [6Jl w here we use calculation methods estab- 
lished in Ref. [23| for forward scattering. At low voltages 
U < 0.17V, quasi-elastic scattering is relevant and we 
take it into account in our numerical calculations by in- 
ner SL valley scattering. This approximations is exact 
for voltages lower than the SL side valley energy gap. 
We can simplify even further the model to an effective 
two- valley model with one central SL valley and one side 
valley by using the approximation of periodic boundary 
conditions for the positions of the potential valleys in 
momentum space. 

In Fig. 2, we show our results for the conductance, the 
differential conductance, and the position and energy av- 
eraged phonon number nr, for certain (m^ma) values 
and lengths L. Here nr mediates the inner SL valley scat- 
tering. Note that an upper limit for 2m!, 2m,2 is given 
by the number of excitable circumverential phonons, i.e., 
2toi,2to2 % d/y/3a where a is the NT interatomic dis- 
tance a w 1.42 A. We obtain a strong increase in the ab- 
solute conductance and differential conductance at high 
voltages (U rj 2V) as a function of mi and m 2 while 
nr is strongly decreasing. The reason for an increase of 
the conductance for larger mi values and fixed 777,2 comes 
mainly from the fact that due to the band edges of the 
side valleys scattering from the central SL valley to the 
SL side-valleys is effectively forward. The backscattering 
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to the central valley is accomplished then by a number 
of different phonons in contrast to the system without 
chiral potential. The growth of the conductance as a 
function of mi is then seen from our discussion in the 
second paragraph which also leads to the explanation of 
the conductance increase as a function of mi . 

Summarizing, we have shown that large chiral unidi- 
rectional superlattice potentials in metallic NTs should 
lead to a large increase of the conductance, the differen- 
tial conductance, and to a decrease of the optical phonon 
temperature at high voltage. We have shown this explicly 



for a chain-like SL potential with symmetric steps. This 
kind of SL potential leads to the best transport results 
over all step potentials. Nevertheless, the main effect 
should be also observed for at least other SL step po- 
tentials with strong chirality. The effect arises from an 
increased number of phonons with different momenta but 
lower electron-phonon scattering probabilities contribut- 
ing to the electron-phonon scattering process. As a result 
of our findings we expect an increase of the applicability 
of carbon NTs as metallic wires. 
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